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Abstract. We analyze numerically how the voltage-current (V — I) characteris- 
tics near the so-called Berezinskii-Kosterlitz-Thouless (BKT) transition of 2D super- 
conductors are affected by a random spatial Gaussian distribution of critical temper- 
ature inhomogeneities with long characteristic lengths (much larger than the in-plane 
superconducting coherence length amplitude). Our simulations allow to quantify the 
broadening around the average BKT transition temperature Tbkt of both the expo- 
nent a in V oc I a and of the resistance V/I. These calculations reveal that strong 
spatial redistributions of the local current will occur around the transition as either / 
or the temperature T are varied. Our results also support that the condition a = 3 
provides a good estimate for the location of the average BKT transition temperature 
Tbkt, and that extrapolating to a — > 1 the ot(T) behaviour well below the transition 
provides a good estimate for the average mean-field critical temperature T c0 . 
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1 Introduction 



As first suggested by Berezinskii [T] and by Kosterlitz and Thouless [2] (BKT) for su- 
perfluids, and later for superconductors by various authors (see, e.g., Refs. |3]j5]), one 
of the main features to be expected in the normal-superconducting transition of a two- 
dimensional (2D) and homogeneous type-II superconductor is the appearance at zero 
applied magnetic field of two critical temperatures: Namely, the transition is split into 
i) the mean-field critical temperature T c o, where it first becomes favorable (in terms 
of free-energy optimization) to form Cooper pairs and vortices, and ii) the so-called 
BKT critical temperature Tbkt < T c o, where it first becomes favorable that vortices 
of opposite fluxoid quantization bind into pairs (vortex-antivortex pairs). Among the 
experimental features revealing this splitting of the transition, probably the most sig- 
nificant is the appearance of a strongly non-Ohmic behaviour in the voltage-current 
V — I characteristics below Tbkt, with the exponent a in V oc I a jumping at Tbkt 
itself to the value a = 3. [5] This feature has been, in fact, commonly used to experi- 
mentally demonstrate the existence of a BKT transition, both in low-T c0 2D structures 
(see, e.g., Refs. (6}jl0]) and in high-T c0 superconducting cuprates (HTSC) (see, e.g., 



Refs. 11-15]). We note that in the case of the HTSC materials even bulk samples are 
expected to undergo a BKT transition, due to the anisotropic 2D-like layered structure 
of these superconductors. [16] In fact, in HTSC with optimal doping the measure- 



ments suggest that the difference T c0 — Tbkt, henceforth noted as Abkt, may be rather 



large, of about 2K. 11 15 This value agrees well with the theoretical predictions, as 
Abkt m &y be approximated by the so-called Levanyuk-Ginzburg criterion |17| that 
also estimates the size of the full-critical (or non-Gaussian) region of superconductivity 
fluctuations above T c q; this full-critical region has been determined on both theoretical 



and experimental grounds to span over about 2K in optimally-doped HTSC. 17 19 



In spite of the fact that most real samples are expected to have some degree of inho- 
mogeneities of the values of T c0 (and of Tbkt), almost no calculations have been done on 
how a spatially-random distribution of T c q- values will affect the BKT non-Ohmic char- 
acteristics. Recently, Benfatto and coworkers |20| proposed a renormalization-group 
study for such situation, but their approach implies to estimate the global resistivity 
of the sample just averaging the ones of the homogeneous domains. This assumption 
could be expected to be adequate only if the current itself is homogeneous in the sam- 
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pie (see also below). In this paper, we use mesh-circuit numerical analyses to study 
the BKT non-Ohmic features that result from considering a 2D type-II superconductor 
having a Gaussian distribution of inhomogeneities of T c0 and Tbkt ; randomly located 
in space and with long characteristic lengths (much larger, in particular, than the in- 
plane superconducting coherence length amplitude). Our analysis allows to obtain the 
evolution of a with temperature, and also shows that significant current redistributions 
occur in the sample as T and I are varied. We believe that our results are applicable 



to measurements (e.g., those in Refs. 11-15]) that show an a(T) jump near Tbkt well 



smoother than predicted by the theory of homogeneous superconductors. 

To get a first glimpse of some of the main difficulties of the proposed problem, let 
us consider in this preface two oversimplified cases of T c -inhomogeneities: i) First, a 
film with domains corresponding to rectangular halves situated with respect to the 
current contacts in series configuration, and ii) the same situation but with domains 
in parallel with respect to the current contacts. In the first case, obviously the current 
passes through both zones without any spatial redistribution as T or / varies, and the 
total resistivity is simply the average of the resistivities of both zones. However, even in 
this uncomplicated case the exponent a of the whole sample will not be just the average 
of the a- values of both zones, as the larger contribution to the total V drop (and hence 
to the global V oc I a behaviour) happens in the zone with larger resistivity. Now we 
consider the situation ii) where the two T c -zones are in parallel configuration. In this 
case, when either T or I are varied, and with them the quotient of the resistivities of 
the two zones, there will be spatial redistributions of the currents. In fact, at some T-I 
combinations, these redistributions will be extreme enough as to become percolating- 
like. Correspondingly, the contribution to a from both zones will be now very different 
to the one in case i). For instance, the global resistivity will greatly differ from the 
average of both zones and the main contribution to a will be given now by the zone with 
lower resistivity. Obviously, a realistic model of a randomly inhomogeneous sample will 
include many domains both in series and in parallel. It will be then nontrivial to know 
if a will be dominated by the zones with higher or lower T c values. Also, a sizeable 
part of the T, / phase diagram will be affected by percolating-like effects. This makes it 
difficult to successfully formulate in a comprehensive T, / range an analytic estimate for 



a in terms of simple averages or even of effective-medium approaches. 21 Because of 



these difficulties, we have chosen numerical simulation methods to analyze the problem. 
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This paper is organized as follows: In section 2 we will briefly summarize the theory 
expressions for the V — I curves of an homogeneous superconducting film near its BKT 
transition. In section 3 we will detail our numerical algorithms and procedures. In sec- 
tion 4 we will present and discuss the resulting V — I curves, and also briefly compare 
them with some previous experimental measurements by other authors, focusing on 



the case of HTSC 12, 13 Section 5 summarizes our conclusions. 



2 Summary of the theoretical approaches for ho- 
mogeneous systems 

2.1 Superconducting contributions to the conductivity near 
the BKT transition at zero applied magnetic field 

To study the V — I characteristics of a 2D type-II superconductor with a random 
spatial Gaussian distribution of Tbkt and T c0 , we will use as starting point the V — I 
expressions for homogeneous superconductors as proposed by Halperin and Nelson 
(HN) in Ref. jH]. To summarize these expressions in a way convenient for our purposes, 
we consider the different T-regimes that appear as T moves from higher to lower values. 

i) Temperatures T > Tbkt- In this T-range the superconducting contribution to 
the electrical conductivity corresponds to the existence of thermal fluctuations of the 
order parameter without vortex-antivortex binding effects, which may be taken into 
account by the Gaussian- Ginzburg-Landau approach (appropriate for T > T c0 + Abkt) 
and the XF-model renormalization-group approach (for Tbkt < T < T c0 + Abkt)- 
HN have proposed [5] an useful interpolation formula that covers the results of both of 
these approaches, with accuracy well sufficient for our present purposes: 

0.37cr n . 2 / 6 A B kt , f , , . 

Act = —r smh W - — (for T > T BK t), (1) 

V 1 ~ BKT 

where Act and a n are, respectively, the superconducting and normal contributions to 
the in-plane electrical conductivity, and b is a dimensionless parameter for which HN 
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do not propose any definite value, stating only that it may be expected to be of the 
order of unity. However, let us already mention here that the value of 60 will be further 
constrained when considering the expressions for T < Tbkt (see point ii below). Note 
also that equation ^ corresponds to an Ohmic Act (if considering an Ohmic a n , as 
will be done in all this work). We also note that we have checked that using the more 
accurate expressions [3}]5|[22] for Act instead of the interpolated equation ([I]) does not 
affect in any significant way the main results presented in this paper. 

ii) Temperatures T < Tbkt- In this T-range, the relevant degrees of freedom for 
the superconducting fluctuations are vortices and antivortices. Again, this region may 
be divided into two, according to the strength of those fluctuations: In the range of 
temperatures closer to T B kt than A B kt the fluctuations are full-critical, while for lower 
temperatures the superconductor follows a conventional Ginzburg-Landau behaviour. 
Following again the ideas of HN we summarize into one common expression the Act 
results for both T-regions at fixed current density j: [6] 

where jo = ek-Q Tbkt / ( ^C«6bkt ) is the Ginzburg-Landau critical current density at 
Tbkt, e the electron charge, k# the Boltzmann constant, h the reduced Planck con- 
stant, d the sample thickness, £ a &BKT = ^afe(0)(T c0 /A B KT) 1 ^ 2 and £ a b(0) the Ginzburg- 
Landau in-plane coherence length extrapolated to respectively Tbkt and T = OK, 
a = max(aT,«j), aj = 3 — l/ln(j/j ) 2 , and ay is: 



3 + W f B A T ~ T (forTBKT-A B KT<T<TBKT), 
V OoAbkT 

(3) 

1 + 2h^f^- (for T < Tbkt - A BKT ). 

^BKT 

Note that «t — > 3 when T — > Tbkt from below, and then at Tbkt it is Act oc j -2 , 
i.e., V oc I 3 (neglecting the small effect of aj and of any other contributions to the 
conductivity). The parameter b\ in equation ^ takes into account the variations of 



the value of Abkt with respect to the purely 2D non-fluctuation GL value. 16 In the 
theoretical simulations done in this work, for simplicity we will neglect such differences 
and take b\ = 1 (except for the comparison with the experimental data shown in figure [6 



where hi has been fine-tuned to the value hi = 2). Note also that to ensure continuity 
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at T = Tbkt — Abkt of equation ^ (and thus of equation (J2J)) it is needed that the 
parameter bo takes the value 60 = 7r 2 /(4&i — 2) 2 , rather than being a somewhat free 
choice as suggested by HN. p\ Finally, note also that the V — I exponent a resulting 
from equation ^ at T < Tbkt — Abkt extrapolates to the Ohmic value at T = T c0 , 
suggesting a possible procedure for the experimental identification of the mean-field 
critical temperature. 

2.2 Other contributions to the conductivity near the BKT 
transition 

To obtain the total conductivity a of the superconductor we must add to the above 
formulas for Aa the contributions from the rest of electrical transport channels in the 
system, mainly the conductivity o n due to the normal-state carriers: 

a = Aa + a n . (4) 

Indeed a n will be negligible against Act for T < T B kt, but for larger temperatures it 
won't be so. Note also that a n may be T-dependent. For instance in the optimally- 
doped HTSC it may be well approximated as inversely proportional to T (as will in 
fact be used in our simulations, see below). We will always consider in this paper that 
a n itself is Ohmic. 

Other contributions to the total conductivity will be neglected in this paper, but 
we note here that in certain specific experimental circumstances they could become 
appreciable: For instance, for T < Tbkt we neglect the Ohmic conductivity that may 
appear at very low intensities when the vortex-antivortex pair breaking processes in- 
volve distances larger than the inhomogeneity size (see, e.g., Ref. [5]). We also neglect 
the non-Ohmic contributions to Aa that may appear above Tbkt due to the super- 
conducting fluctuations not related to the vortex-antivortex correlations. The latter 
contributions have been thoroughly explored previously (see, e.g., Refs. 
today are well known to be significative only at electrical fields much larger than those 
considered in our study (see, e.g., figure (1) of Ref. [27]). Finally, we also mention that 
we neglect the indirect contributions to the paraconductivity above the transition, such 
as the so-called Maki-Thomson and density-of-states contributions. [19] It is today quite 



23-26n and 
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well accepted that such contributions are negligible in the case of HTSC 19,22,28 al- 



though the situation is not as clear in the case of low-T c superconducting thin films. 29 



We have checked however that including any of such contributions in our simulations 
does not qualitatively affect the main results presented in this paper. 



3 Procedure for the numerical simulation of the 
V — I characteristics near the BKT transition of 
an inhomogeneous superconductor 

Our aim is to obtain the V — I characteristics of a 2D superconductor composed by 
randomly-located domains, each domain having its own single T c q and Tbkt and fol- 
lowing the V — I characteristics described in the previous section. For this purpose, in 
the spirit of the finite-element methods we model the inhomogeneous superconductor 
as a N x N square mesh of resistors (see figure [I] and figure [2]). We randomly assign to 
each node of the mesh a different T c0 , and to each resistor the T c0 of its corresponding 
left-lower node. The T c0 distribution is Gaussian with mean-value T c0 and full-width 
at half-maximum AT c0 . The difference T c0 — Tbkt is held constant for all resistors, so 
that the distribution of Tbkt follows the one of T c o and is Gaussian with mean- value 
Tbkt and full- width at half-maximum also AT c o. We include also in our model an 
external circuit composed by a current source connected to opposite borders of the 
sample with zero-resistance contacts. When referring to the results of our simulations, 
by / we mean this external bias current and by V the voltage drop between those 
opposite contacts. The sample is considered to have width and length w and thickness 
d. 

Because the resistivity of each resistor depends on the local current passing through 
it (see previous section), the mesh equations that result from that modeling are nonlin- 
ear and in general do not admit analytic solution. They have to be solved using numer- 
ical methods, of which we use Newton-type iterations. For those iterations to succeed, 
it is crucial to start them from initial values not too far from the solution. Thus, to 
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calculate each V — I curve at a fixed temperature, we applied the following algorithm, 
that proved itself to be well adapted to the non-Ohmic features of the BKT transition 
(as it solves the instability problems that we found trying other simulation strategies): 
Our analysis starts by considering first a high temperature T start 3> T c0 (where the 




Figure V. Some examples of the evolution of the current distributions obtained in our 
simulations when the temperature is varied and the bias intensity is kept constant. 
The simulation sample is composed by a 10 x 10 mesh of resistors corresponding to a 
random spatial Gaussian distribution of mean-field critical temperatures with average 
T c o and full-width at half-maximum AT c o, with the value of Abkt = T c q — Tbkt being 
the same for all resistors. In the pictured example, T c0 = 100K, Abkt = 2K (so that 
Tbkt — 98K), AT c0 = 2K, the sample width and length is w = 10 _3 m and its thickness 
is d = lOOnm. We also used (see main text) b± — 1, a a — 10 S T~ 1 K/Qm and £ a b(0) = 
lnm. The bias current / is the one provided by the external current source, which 
is connected to opposite borders of the superconductor with zero-resistance contacts. 
Note the strong spatial redistribution of currents inside the superconductor as T is 
varied through the transition, indicating that at each temperature different regions 
of the superconductor determine the global resistivity, and its Ohmic or non-Ohmic 
character. These redistributions are specially intense near the average BKT transition 
temperature where the current flows through only a few less resistive paths. It may 
be also observed that for T ^> T c o almost all of the current flows longitudinal and 
uniformly, while for T < Tbkt transversal current paths are significant. 
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system is Ohmic and easy to solve) and a high bias intensity I s tart = {wd/10)jo- The 
analysis then evolves keeping the bias intensity constant but lowering the temperature, 
evaluating at each T-step the mesh equations by means of Newton iterations with 
starting point the final result of the previous step. The T-decrement separating each 
step is adaptively updated during the simulation, so that voltages are not allowed to 
vary beyond 0.1% between steps. Once reached the target temperature at which the 
V — I curve is to be calculated, the temperature is fixed and then the bias intensity 
/ is varied, again adaptively, iteratively solving the mesh equations at each step and 
storing the results. Each J-step uses as starting values for the Newton iterations the 
results obtained in the previous calculated step. As an additional measure to avoid 
instabilities in the convergence of the Newton iterations, when needed our program 
smoothes over 0.05K the V(T) behaviour of the individual resistors at their local Tbkt 
temperature, making its V(T) evolution continuous (but still very rapidly varying). 
This T-widening of the BKT transition is negligible in any case against the one due 



T = 98K = Tbkt 
AT c0 = 2K 

1 = 0.1 mA 



1/10 



1/100 




1= 1 mA 




I = 10 mA 




Figure 2: Some examples of the current redistributions obtained in our simulations, 
when the bias intensity is varied and the temperature is constant. The simulation 
sample is the same as in figure [T} The fact that different regions of the inhomogeneous 
sample contribute to the electric transport as / changes indicates that the log-log slope 
of V(I) may vary as / is varied, in spite of the power-law-like V oc I a behaviour of the 
homogeneous case. This change in slope is appreciable in the figures [31(a) and|3| 
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to the inhomogeneities considered in our simulations and we checked that doubling or 
halving it does not change our final results in any appreciable way. 

All of the results presented in this paper correspond to 10 x 10 cell meshes. We 
have checked that runs of the simulation using different random distributions but with 
the same statistical parameters (T c0 , Tbkt and AT c0 ) produce similar outcomes. Com- 
putation time for a single V — I curve at fixed T is of about three days for a 10 x 10 
mesh in current desktop computers. 



4 Results of the simulation 



Let us now comment on the results obtained when applying the method described in 
the previous section to compute the V — I curves around the BKT transition, using 
parameter values typical of HTSC film samples. In particular, we have used T c0 = 
100K, Abkt = 2K (therefore T BK t = 98K), a n = K^T^K/fim, £„ 6 (0) = lnm, d = 
100 nm and w = 10 _3 m, and we have constructed simulated samples from AT c q = 2K 
to 4 K. 

In figure [T] and figure [2] we show some examples of the current distributions within 
the mesh circuit, obtained for representative values of T and /. The most important 
feature observed in these results is that significant spatial redistributions of the currents 
may occur when either T or / are varied. This is a consequence of the fact that the 
resistance of each mesh element relative to the resistance of the other elements will be 
dependent on T and /, and therefore the current paths will be also dependent on both 
of these variables. Due to these redistributions, changes of T and / will also change 
what are the portions of the sample that dominate the global voltage drop, and its 
Ohmic or non-Ohmic character. 

We may also note that, as in fact it was already commented in the Introduction, 
the behaviour of the system in general will be intermediate between the simplest cases 
of considering all resistors in series (in which case the global resistivity is given by 
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the average of the resistivity of all the elements) or in parallel (in which case the 
minimum resistance will dominate the transport properties), being the proximity to 
each situation dependent on / and T. Interestingly, for temperatures above T c0 the 
current distribution becomes especially simple: As it can be seen in figure[I]and figure [2j 
at T > T c o almost all of the current flows longitudinal and uniformly. So, at those 
temperatures the global resistance may be approximated by the one of a row, which in 
turn corresponds to the average resistance of the resistors on it. As the temperature is 
lowered we see however that the current path geometry is no longer that simple, nor 
constant with T or /. Even for T < Tbkt the current does not flow longitudinally, and 
transversal current paths remain significant. 




10-4 -1 0-3 10-2 -1 0-4 10-3 10-2 

I (A) I (A) 



Figure 3: V — I results at various constant temperatures around the BKT transition 
as obtained in our simulations in two planar superconductors with (a) AT c0 = 2K and 
(b) AT c o = 4K. In both cases we have used parameter values typical of HTSC films, 
namely T c q = 100K, Tbkt = 98K, Abkt = 2K, thickness lOOnm, width and length 
10 _3 m, a n = 10 8 T _1 K/fhn, and £ a &(0) = lnm. We have also used b\ — 1 (see main 
text). The log-log slope of these V — I results corresponds to the exponent a (see also 
figure [3]). The T = T B kt isotherm (thickest line) presents a = 3, while the T = T c0 
isotherm (upper line) is Ohmic (a = 1). Note that non-Ohmic behaviour appears in 
these inhomogeneous superconductors already above Tbkt- Note also that the log-log 
slope may vary with 7 or V. In this work the results reported for a(T) correspond to 
the Grange 3 x 1(T 6 to 3 x 1(T 5 volt. 
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In figures[3](a) and[3[b) we show the V — I curves that result from these simulations. 
As mentioned previously, here V and I correspond to the global values, i.e., those in 
the external bias circuit. The log-log slope of these curves corresponds to the exponent 
a in the V oc I a dependence. It is evident in figures [3]^a) and ^b) that a depends on 
temperature. As it could be expected, Ohmic behaviour (i.e., slope unity) is obtained 
for temperatures well above Tbkt, while much larger log- log slope is obtained well 
below that temperature. For temperatures close to Tbkt the change in behaviour is 
not discontinuous, being instead somewhat broadened by inhomogeneities, although the 
T-range where that broadening occurs is significantly smaller than AT c q. This change 
in the log- log slope may be seen more accurately in figure |4j where a(T) is plotted for 
all of the simulated samples, together with the theoretical a(T) corresponding to the 
homogeneous case. To obtain this figure |4j a was calculated trough a power-law fit to 
our V — I results at voltages around 10~ 5 volt (in particular 3 x lCT 6 volt < V < 3 x 




Figure 4: The exponent a in V oc I a extracted as the log-log slope of the simulation 
results shown in figures |3^a) and[3|b), in the range 3 x 10~ 6 to 3 x lCT 5 volt. The 
homogeneous case, AT c q = OK, is also shown for comparison (see section 2). Note that 
the condition a = 3 provides a good estimate for T B kt, and that extrapolating the 
ce(T) behaviour for T < Tbkt — Abkt (=96K) to the Ohmic value, a = 1, provides a 
good estimate for T c0 . 
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lCT 5 volt). The reason why it is necessary to specify a voltage range for the obtainment 
of a is that, as may be noticed in figures |3|a) and[3|b), in the inhomogeneous samples 
for temperatures close to T B kt the V — I dependence is not perfectly power-like, but 
rather the log-log slope depends also on the applied current. The cause behind this 
fact is the existence of local current redistributions, shown in figure [I] and figure [2] (i.e., 
changes in / vary the region of the sample where the voltage drops, and thus its global 
Ohmic or non-Ohmic character). 

Other important feature to be observed in figure [4] is that two of the most common 
criteria used by experimentalists [6 15 to determine from the ot(T) plots the BKT and 



mean-field critical temperatures remain essentially valid in spite of the inhomogeneities, 
if we apply those criteria to determine now their average values, Tbkt and T c0 : In 
particular, the condition a = 3 provides a good estimate for Tbkt, underestimating it 
only slightly (in particular the deviation from the exact value is well smaller than AT c o). 
Secondly, the a(T) dependence for temperatures T < T B kt — Abkt extrapolates to the 
value a = 1 at T = T c q with excellent accuracy, indicating then that such extrapolation 
procedure remains adequate in real inhomogeneous samples to determine T c0 . 

The effects of the inhomogeneities in the dc electrical transport properties become 
more apparent in the R(T) curves obtained at different fixed external currents I. Here 
we define R as simply V/I. In figures [5][a) and|5^b) we show the R(T) results obtained 
from our simulations at different fixed I values, for samples with different AT c q values 
and using a logarithm axis for the resistance. It may be seen in these figures that the 
application of a finite current broadens the logi?(T) tail in the lower part of the super- 
conducting transition, and that this happens to a larger extent as AT c q is increased. It 
is also easily noticeable that the non-Ohmic behaviour sets in at temperatures above 
Tbkt- We conclude that the inhomogeneities are detectable over a larger T-range on 
the amplitude of the resistance than on the a exponent. 

Let us now discuss how the results obtained with our simulation procedure compare 
with some of the experimental data obtained by earlier authors measuring the V — I 
characteristics near the BKT transition. Specifically, we will use for these comparisons 



the a(T) measurements of Refs. 12 and 13 , as in these works a{T) was extracted 



using a voltage criterion similar to the one used in our simulations. These measurements 
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were performed in Ref. [12] in T^B^CaC^Og 700nm-thick films (see in particular 
figure 3(a) of that work) and in Ref. [13] in YB^CusOy-^ 120nm-thick films (see figure 2 
of that work). We show in our figure [6] a comparison between these experimental a(T) 
and our simulation results. To be able to gather together in a single representation 
the two samples in spite of their different critical temperatures, we have chosen for 
the horizontal axis the normalized quantity (T — Tbkt)/Abkt- We employed for each 
sample the T B kt that results from applying the condition «(T B kt) = 3, and the T c0 
that results from extrapolating to a — > 1 the low-temperature a(T) data. We obtained 
T BKT = 99.0K and T c0 = 100.2K for Ref. |l§ and T BKT = 83.45K and T c0 = 85.95K 
for Ref. 13 . In our simulations we have used the same parameter values as for the 
simulations shown in figures 1 to 5, except for AT c q = 4K and b± = 2. These latter 
values were found to be the ones producing a better agreement between those data and 



i i i i ri n i i i i r 




T (K) T (K) 



Figure 5: Resistance versus temperature at various constant bias currents as resulting 
from the quotient V/ 1 in the same simulation runs as in figures |3](a) and [3^b) for 
superconductors with (a) AT c0 = 2K and (b) AT c0 = 4K. The logarithmic scale allows 
to better appreciate that the tail of the resistive transition is broadened when the 
current is increased and that the effect is larger for the samples with a larger AT c o. 
Note also that non-Ohmic behaviour sets in already above Tbkt- 
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our simulations (note that in the case of b\ its value is governed mainly by the data for 
T < Tbkt — Abkt, outside of the region significantly affected by inhomogeneities) . As 
may be seen in the figure |6j the agreement between experiments and simulation is rather 
satisfactory, despite the necessary crudeness of some of our approximations (perfectly 
Gaussian distribution of the T c o-inhomogeneities, uniform Abkt value, inexactness of 
any finite-element method, etc.). 




(T-TbktVAbkt 



Figure 6: Comparison between the exponent a obtained experimentally in Refs. 12 
and [13) in HTSC (circles and squares) and our simulations (solid line). The ex- 
perimental data were obtained by means of V — I measurements in Tl 2 Ba 2 CaCu208 
700nm-thick films 12 and in YBa2Cu3Q7_g 120nm-thick films Jl3|. In this comparison 
we employed for T B kt the values that result from the condition a = 3, and for T c0 
the values that result from extrapolating to a — > 1 the low-temperature ct{T) data. 
We obtained T BKT = 99.0K and T c0 = 100.2K for Ref. |l| and T BKT = 83.45K and 
T c0 = 85.95K for Ref. pi. We also used the values AT c0 = 4K and 61 = 2, that 



produce the best agreement with these a(T) data. 
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5 Conclusions 



We have analyzed numerically the effects of a random spatial Gaussian distribution of 
critical temperature inhomogeneities with long characteristic lengths on the voltage- 
current V — I characteristics of a type-II planar superconductor near the Berezinskii- 
Kosterlitz-Thouless (BKT) transition. The simulations allow to quantify the broad- 
ening around the average BKT transition temperature Tbkt of both the exponent a 
in V oc I a and of the resistance V/I. These calculations reveal that strong spatial 
redistributions of the local current will occur around the transition as either I or T are 
varied. Our results also support that the condition a = 3 provides a good estimate for 
the location of the average BKT transition temperature Tbkt, arid that extrapolating 
to a — > 1 the ot(T) behaviour well below the transition provides a good estimate for the 
average mean- field critical temperature T c0 . These results are in good agreement with 
some experimental measurements of the exponent ot(T) obtained by earlier authors on 



HTSC films. 12,13 
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